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Abstract 



We define a block observable for the g-state Potts model which exhibits an 
intermittent behaviour at the critical point. We express the intermittency 
indices of the normalised moments in terms of the magnetic critical expo- 
nent (3/v of the model. We confirm this relation by a numerical similation 
of the q = 2 (Ising) and q = 3 two-dimensional Potts model. 
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1 Introduction 



Recently, the concept of intermittency has been intoduced in the field of 
equilibrium critical phenomena as a tool for studying local fluctuations of 
a system undergoing a second order phase transition Q ||]. For example, 
in the Ising model, as the the magnetisation fluctuates without limit at 
the critical point, one expects that the normalised moments of this observ- 
able, measured on sub-blocks of the system, scale with the size of the block, 
with exponents that may depend on the order of the moment, the so-called 
intermittency indices. On the basis of real space renormalisation group ar- 
guments, Satz || established that, for the Ising model, the intermittency 
indices are directly connected to the magnetic critical exponent. However, 
the normalised moments of the magnetisation used in Satz's argument are 
not suitable for a numerical simulation, as the net magnetisation, which 
enters in the denominator of the moments, is almost zero on a (large) fi- 
nite lattice. In subsequent numerical simulations, Bambah et. al.||j and 
Gupta et. al.j5| ; using a new observable which respects the Z2 symmetry 
and remains finite on finite lattices, present clear evidence for an intermit- 
tent behaviour in the two-dimensional Ising model. However, they do not 
observe the expected relation between the intermittency exponents and the 
critical ones. Burda et. al. @ suggest that this is a consequence of their 
choice for the fluctuating observables, among which the moments have only 
a subdominant scaling behaviour in the critical regime. 

Therefore, the problem of finding convenient block observables which 
have well-defined moments on finite lattices and which exhibit the expected 
scaling law, remains open. In this paper we address this question in the 
more general framework of the g-state Potts model. 

In section 2, we define the fluctuating observable for which we derive 
the intermittency indices in terms of the magnetic exponent by generalising 
the Satz's argument §. In section 3 we give the details of the numerical 
simulation and discuss the results. Conclusions are drawn in the last section. 

2 Definition of the observables 

The most convenient definition of the moments for studying critical spin 
systems has been thoroughly discussed in ref [jj], |5|. Following the authors 
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of ref. P| we use the so-called standard block moments of order p defined as 

M < k a >p 



where L is the lattice size, £ the cell size, M = (L/£) d the number of cells 
and k a the block observable defined on the a th cell. The symbols <> stands 
for the thermodynamical average which will be taken at the bulk critical 
temperature T c defined by K c = J/kbT c = ln[l + ^/q] for the g-state Potts 
model. We recall that the intermittency indices X p are defined according to 
the behaviour of the moments with respect to the block size £ : 

F p {L;£)~(j-) for-^0 
Let us now define the block observable k a . 

We consider a system of Potts spins Sj = 0, 1, . . . , q— 1, interacting with 
their nearest neighbours via the hamiltonian 

-(3H = KJ2 *s iSj 

<i,j> 

on a two-dimensional square lattice of linear size L. For a particular con- 
figuration, we denote by Q the spin value which occurs the most frequently 
and to each site i we assign the variable n« = S Si Q. We now define 

_ qrii - 1 



q-l 
which is such that 

Jd < ki >=< kl >=< m i > 

i 

is the order parameter for the Potts model [|7j. For the Ising model (q = 2) 
this definition reduces to 

ki = sgn(m/)<7i 

where mj is the Ising magnetisation of the configuration of Ising spins {o~i = 
±1}. 

We then define the block observable to be 

1 qrig ' — 1 

ka = gd ki = q _i ( 2 ) 
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where is the fraction of spins in the cell a having the value Q || . 

Unlike the original definition of ref pi, on a finite lattice, < k a > is 
non zero at the critical point, (< k a > = < \mj\ > for the Ising model) 
allowing the moments of eq. (2.1) to be measured in a numerical simulation. 
Moreover, one can extend to this observable the renormalisation argument 
which predicts the intermittent behaviour of the moments. 
Let us associate to the blocking procedure the following renormalisation 
prescription : we define a block spin s a according to the majority rule in the 
a th cell, associate the corresponding renormalised site variabld^ Uq = 5 Sa Q 
and then the renormalised quantity k a 

_qn^-l 

fun- — 



q-1 

Since k a is related to the magnetic scaling operator, we expect the following 
renormalisation relation 

where Z{€) is the renormalisation factor depending only on the block size t. 
Therefore 

<!>■>■> [Z (()]"< I 9r ^ ) ~ 1 ) > (3) 



(4) 




[q-lY / (q-l)P 



where we have used {uq'Y = . Reintroducing k a in the left-hand side 
of eq.(2.4), we get the relation 

< kl >= [Z{l)f {A p <k a > +B p ) (5) 

where the coefficients A p and B p are given exactly by 

a i (- l ) (p+1) \ r «x 



1 For a given configuration, the value of Q is preserved by the blocking procedure 
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According to the standard renormalisation analysis, we expect Z{£) ~ 
£~^(d-2+r]) _ £-f3/u _ Moreover, < k a > and < k a > are the order parameter 
of the system of size L and of the rescaled system of size L/£, respectively. 
At the bulk critical temperature, these quantities behave according to the 
finite size scaling law 



<k a > ~ L-M" 
<k a > ~ (L/£)~^ u 



(8) 
(9) 



Thus, for £ S> 1 (in lattice units) in order to sum up enough degrees of 
freedom for the renormalisation arguments to be valid and £ <C L to avoid 
finite size effects, we expect the moments F p to behave like 



Mi 



+ B' 



(10) 



where the coefhcients A' and 5' are proportional to the coefficients A p and 



B p defined above. 



For the Ising model (q = 2), eqs (2.6) and (2.7) give A p = for p even 
and Bp = for j> odd, leading to the scaling law, derived previously by 
Satz H 



F p (L;£) 



n-(p-i)l 



for p even 
for p odd 



(11) 
(12) 



For the three-state Potts model such a cancellation does not occur and 
the leading behaviour (£/ 'L)~ p P/ u is affected by a corrective term as shown 
in eq. (2.10). However, similar corrections also appear in the scaling laws of 
eqs (2.9), which, when taken into account, modify the exact exponent (3 /is 
in eq.(2.10). Therefore, for the three-state Potts model, we expect a more 
general form 



F P (L- 



L 



+ b; 



(13) 



where Ap, B' p ' and w are unknown parameters. 
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From this behaviour we deduce the intermittency indices 



For the Ising model 




1) 



for p odd 



for p even 



(14) 



For the three-state Potts model 



^p 
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3 Numerical tests 

We have performed a simulation of the q = 2 and q = 3 Potts models at 
the bulk critical temperature, using the Swendsen Wang dynamics, on two- 
dimensional lattices of size up to 256 x 256. For the largest lattice sizes, the 
data were taken on four independant runs of 2 x 10 4 Monte-Carlo lattice 
updates, with a spacing between two consecutive measurements of five MC 
lattice updates. The error analysis is realised on this total sample of 1.6 x 10 4 
independent measurements, leading to quite small error bars (of the order 
of the width of the data points). In order to check the reliability of our 
simulation, we have tested the scaling law for the order parameter, < m q > 
(eq. 2.8), for lattice sizes ranging from L = 16 up to L = 256. The results 
are displayed in figure 1, with a linear fit in log scale, giving (3/v = 0.124(5) 
for the Ising model (the exact value is 1/8=0.125) and (3/v = 0.131(7) for 
the three-state Potts model (the exact value is 2/15=0.133). 

We then determine the intermittency indices X p . We measure the mo- 
ments for p = 2 to p = 5 and for £ = 2, 4, 8, ... , L. Notice that the values of 
the moments for £ = L approach, as L —* oo, the ratio of universal critical 
amplitudes. For instance, the ratio [F2(L, £ = L)] 2 /F<i(L, £ = L) has been 
thoroughly studied || and found to have the value 0.85622 in the Ising 
model |K|. Our result is 0.856(4) providing a further check of our data. 

In figure 2 we show the moments for the Ising model on a 256 x 256 
periodic lattice. The behaviour of F p depending on the parity of p depicted 
by eqs (2.11-12) is quite visible on the data. According to these equations, 
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we have performed a power law fit of the moments for 4 < I < L/2, shown 
as the straight lines in figure 2. We observe that the odd moments are 
remarkably well fitted by the power law, whereas the even moments exhibit 
deviations from this leading behaviour. Therefore, we proceed to a double 
determination of the exponents : first, from the exact power law of eqs.(2.11- 
12) and second, from a corrected parametrisation of the form of eq.(2.13). 
The results of both fitting methods for the largest lattice size L = 256 are 
shown in the last two lines of table I. The difference between the values 
corresponding to the same moment give an estimate of the systematic error 
on the exponents. In addition, we give the result of the power law fit for the 
L = 128 lattice, which shows the stability of this determination. 





P 


L 


2 


3 


4 


5 


128 


0.111 


0.124 


0.117 


0.123 


256 


0.114 


0.124 


0.119 


0.123 




0.135 


0.126 


0.129 


0.122 



Table I The exponents X p /p (p even) and X p / (p — 1) (p odd) for the Ising 
model on lattices of size L = 128, 256 and for the moments of order 
p = 2, 3, 4, 5. The exponents given in the bottom line correspond to a 
fit of the L = 256 moments with a corrected scaling law. The expected 
exact value is (3/v = 0.125. 



The agreement of these results with the prediction of eqs.(2.14) j3/v = 
0.125, is qualitatively good for the even moments and excellent for the odd 
ones. 

We have repeated the same analysis for the three-state Potts model. The 
moments are shown in figure 3 for the largest lattice size L = 256. The fit 
corresponds to the corrected power law fit of eq.(2.13) 



F p (L;£) = 




where x, w, a, b are free parameters. The results are shown in table II for 
the lattice size L = 256. 
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L 


P 


2 


3 


4 


5 




X 


0.141 


0.124 


0.131 


0.125 






± 0.015 


± 0.009 


db 0.009 


db 0.011 


256 


w 


0.359 


0.270 


0.307 


0.235 




a 


0.508 


0.581 


0.652 


0.772 




b 


0.516 


0.503 


0.521 


0.560 


128 


X 


0.144 


0.119 


0.128 


0.124 






± 0.022 


db 0.017 


db 0.020 


db 0.011 



TABLE II 



Table II The exponents x = X p /p for the moments of order p = 2, 3, 4, 5 for the 
three-state Potts model on the L = 256 lattice. The expected exact 
value is (3jv = 0.133. We show the fitted values of the other variables 
entering the parametrisation of the moments F p (z) = z~ px {az w + b) 
where z = i/L. x values from the L = 128 data are reported at the 
bottom line. 

The errors quoted are the statistical ones corresponding to 95% con- 
fidence level. They are quite large due to the freedom allowed by our 
parametrisation. Actually, we observe that, if we set the exponents to the 
exact value, x = 0.133, the minimum x 2 does not change significantly with 
respect to its best-fit value. We give the other variables of our parametri- 
sation in Table II. They are only weakly dependent on the order of the 
moment, but the exponents w differ significantly from x, indicating, as ex- 
pected, a strong contamination of the behaviour of eq.(2.10) by corrections 
to the scaling law of eq. (2.9). The exponents obtained from the fit of the 
L = 128 data are reported at the bottom line of table II, showing again the 
stability of their determination. 

Nevertheless, the agreement with the expected intemittency indices is 
again quite good. 

4 Conclusions 

We have defined a block observable for the g-state Potts model which ex- 
hibits an intermittent behaviour at the critical point of the model. The 
moments of this observable are measurable on finite lattices and have well- 
defined scaling behaviour in terms of the block-size. By a numerical simula- 
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tion, we confirm the previously observed 0, |B[ intermittent behaviour, and 
we find good agreement between the measured exponents and the magnetic 
exponent (3/v, as predicted by renormalisation group arguments. 

Although a critical spin system is presumably not multifractal, a com- 
plete description of its scaling properties requires several exponents. Besides 
the thermal and magnetic critical indices which drive the fluctuations of en- 
ergy and of magnetisation, the geometrical aspects of the critical clusters 
of spins are characterised by other exponents, not related to the thermo- 
dynamics ones [0, m, U. For this reason, the intermittency phenomenon 



in a critical system may assume several aspects depending on the definition 
of the block observable. This may explain the results of ref ||, where the 
measured intermittency indices coincide with the fractal dimension of the 
clusters of spins. However, further theoretical and numerical studies are 
needed in order to investigate the geometrical origin of intermittency. 
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Figure Captions 



Figure 1 The order parameter for the q = 2 (Ising) and q = 3 Potts model, 
measured at the bulk critical temperature as a function of the lattice 
size. The lines result from a power fit, giving the exponent (3/v = 
0.124(5) for the Ising model (exact value 0.125) and (3/v = 0.131(7) 
for the three-state Potts model (exact value 0.133). 

Figure 2 The moments of order p = 2, 3, 4, 5 for the Ising model as a function 
of the block size i. The straight lines result from a power law fit. 

Figure 3 The moments of order p = 2, 3, 4, 5 for the three-state Potts model 
as a function of the block size I. The curves result from a corrected 
power law fit. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9303097vl 



This figure "fig2-l.png" is available in "png" format from: 
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